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摘要 能量 扫描 理论 合作 组 (Beam Energy Scan Theory Collaboration ,BEST) 的 目标 是 建立 一 个 可 以 描述 美 医 
相对 论 重 离子 对 撞 机 (Relativistic Heavy-Ion Collider, RHIC) 上 第 二 期 能 量 扫描 实验 的 动力 学 框架 。 该 实验 有 
可 能 找到 强 相 互 作用 相 图 上 的 临界 点 。 本 文 总 结 自 2016 年 起 ,BEST 合 作 组 取得 重要 进展 ,并 对 未 来 探索 
密度 区 相 图 做 了 展望 。 

关键 词 ”量子 色 动力 学 (Quantum Chromodynamics,QCD) 相 图 , 相对 论 重 离子 对 撞 , QCD 临界 点 , BRIE, 动力 
学 模型 
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The BEST framework for exploring the QCD phase diagram: progress summary 


YIN Yi 


(Institution of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China) 


Abstract Exploring the quantum chromodynamics (QCD) phase diagram at finite bayron density regime through 
the beam energy scan (BES) program at the relativistic heavy-ion collider (RHIC) is one of the key frontiers in high 
energy nuclear physics. The high precision data anticipated from the second phase of the BES program would 
potentially enable the discovery of the conjectured QCD critical point, a landmark point on the phase diagram. In this 
paper, the progress made by the beam energy scan theory (BEST) collaboration, which was formed with the goal of 


providing a theoretical framework for analyzing data from BESII, is reviewed. In addition, the challenge of 
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investigating the QCD phase diagram with future facilities is discussed. 


Key words QCD phase diagram, Relativistic heavy-ion collisions, QCD critical point, Fluctuations, Dynamical 
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QCD) 是 描述 夸克 胶 子 相互 作用 的 基本 理论 。 


a 色 动 力学 (Quantum Chromodynamics， 


不 断 取得 习 
研究 WAR, ES 
极端 高 温 高 密 物 质 (QCD 物质 ) 属 性 是 量子 色 动 力 QGP 是 平滑 
学 研究 的 重要 前 沿 。 过 去 20 年 ,美国 相对 论 重 离子 


要 进展 ,例如 最 新 第 


154~158.4 MeV。 这 日 


性 原理 格 点 QCD 


E 子 化 学 势 为 零 时 ,从 强 子 气体 到 
骨 相 变 (Crossover transition) , 相 变 温度 在 
描述 的 都 是 习 重 EE 子 密度 较 低 时 


对 撞 机 (Relativistic Heavy-Ion Collider, RHIC) 以 及 QCD 物质 的 属性 。 学 科 当 前 的 重要 热点 是 系统 研 
欧洲 大 型 强 子 对 撞 机 (Large Hadron Collider, LHC) F QCD 物质 在 中 高 密度 区 的 属性 六 。 


上 的 相对 论 重 离子 对 撞 实验 ,揭示 了 一 系列 QCD 物 


质 的 新 颖 性 质 。 比 如 ,处 在 解禁 闭 状态 下 的 夸 
子 等 离子 (Quark-gluon Plasma,QGP) 被 发 现 是 一 


流动 性 极 好 的 近 理 想 流 体 。 与 此 同时 ,理论 太 


TE 的 临界 点 CCritical point) 以 及 其 
白 点 表示 ,红色 点 对 应 不 同 对 撞 能 量 下 火球 


究 也 freezeout 大 致 在 相 图 
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图 1 是 我 们 现在 对 QCD 相 图 的 认 知 。 猜 想 中 


可 能 存在 的 位 置 由 


图 中 还 展示 了 第 二 
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期 能 量 扫描 实验 (Beam Energy Scan,BES-ID 履 盖 的 
相 图 区 域 。 可 以 看 到 ,临界 点 是 QCD 相 图 的 地 标 。 
很 多 QCD 模型 计算 都 表明 , 随 着 重子 密度 的 增加 ， 
相 边 界 将 从 平滑 相 变 过 渡 到 不 连续 的 一 级 相 变 。 与 
之 对 应 ,一 级 相 变 在 相 图 上 终止 的 点 被 称 作 QCD 临 
界 点 ,但 是 临界 点 是 否 存在 ,现在 还 无 法 通过 第 一 性 
原理 计算 。 因 此 ,寻找 临界 点 以 及 一 级 相 变 的 信号 ， 
是 现 阶段 重 离子 对 撞 实 验 的 重要 目标 。 

2019~2021 年 ,RHIC 进 行 了 BES-I 实 验 , 对 撞 
模式 的 能 量 点 覆盖 19.6~3 GeV ,大 致 对 应 重子 化 学 
势 200~400 MeV。 在 这 样 一 个 能 量 区 间 ,一 系列 对 
临界 点 敏感 的 观测 量 在 第 一 期 能 量 扫 描 实 验 中 呈现 
出 有 趣 的 表现 。 由 于 BES-I 实 验 数据 的 极 高 精度 ， 
该 实验 有 潜力 将 以 前 实验 得 到 关于 临界 点 的 线索 ， 
变 成 其 存在 的 科学 证 据 。 但 是 , 重 离子 对 撞 实 验 产 
生 的 火球 有 着 复杂 的 动力 学 演化 ,从 实验 结果 中 提 
取 临 界 点 的 信息 ,需要 发 展 相应 的 理论 框架 。 在 这 
样 一 个 背景 下 ,由 美国 能 源 部 资助 ,能 量 扫描 理论 合 
作 组 (Beam Energy Scan Theory Collaboration, BEST 
合作 组 ) 于 2016 年 成 立 ,资助 期 为 5 年 。BEST 合 作 
组 的 主要 目标 是 建立 一 个 定量 的 动力 学 框架 ,用 于 
描述 第 二 期 BES 实验 的 数据 ,回答 在 BES 能 区 是 否 
处 在 临界 点 的 问题 。 除 此 之 外 , 这 个 动力 学 框架 也 
可 以 应 用 于 提取 有 限 密度 区 QCD 物质 其 他 重要 
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图 1 QCD 相 图 示意 外 
Fig.1 A sketch of the QCD phase diagram”! 

2016~2020 年 , BEST 合作 组 在 构建 前 述 动力 学 
框架 方面 取得 了 一 系列 重要 进展 。2020 年 资助 期 
结束 之 后 ,合作 组 成 员 之 间 的 相互 合作 也 十 分 活跃 ， 
见 文献 [2] 的 总 结 。 本 文 将 综述 BEST 合作 组 在 
QCD 相 图 研究 取得 的 进展 ,以 期 对 未 来 研究 提供 有 
益 的 参考 ,参看 相关 综述 文献 [3-4]。BEST 合作 组 
的 另 一 个 目标 是 定量 研究 手 征 磁 效 应 (CChiral 
Magentic Effect,CME) 。 对 这 方面 研究 活动 的 评述 
不 在 本 文 的 讨论 范围 之 内 ,有 兴趣 的 读者 可 以 参考 


趣 成 果 , 可 参看 文献 [8-9]。 
对 于 LHC 能 区 以 及 RHIC 最 高 对 撞 能 量 ,在 
2016 年 之 前 已 经 有 一 套 基 于 相对 论 流体 力学 的 成 
熟 动力 学 框架 来 描述 对 撞 产 生 的 火球 的 演化 。 该 杠 
架 用 尽 可 能 少 的 参数 描述 了 一 系列 实验 结果 。 把 该 
框架 拓展 到 BES 能 区 ,需要 解决 的 主要 问题 包括 : 

1) 初 始 条 件 : 随 着 核 核 碰撞 能 量 的 降低 ,其 初始 
状态 会 发 生 一 系列 定性 的 区 别 ,比如 夸克 自由 度 占 
的 比重 将 增 大 ,而 对 于 高 能 量 碰撞 ,自由 度 主 要 是 胶 
子 。 在 这 些 磁 撞 中 ,保留 在 火球 中 的 重子 数 也 在 增 
加 ,其 分 布 对 重子 stopping 的 机 理 十 分 敏感 。 因 此 
需要 构建 考虑 了 这 些 重要 因素 的 初始 条 件 模型 。 

2) 状 态 方程 :由 于 符号 问题 , 格 点 计算 无 法 提供 
重子 密度 较 大 时 的 状态 方程 。 更 为 重要 的 是 ,如 果 
临界 点 被 BES 能 区 覆盖 ,作为 定量 模型 输入 的 状态 
方程 必须 包含 临界 点 的 信息 。 

3) 流 体力 学 :需要 演化 重子 密度 ,奇异 
(strangeness ) fii #4 RE“ H EH RE 

4) 临 界 涨 落 的 演化 :临界 涨 落 增强 是 临界 点 存 
在 的 重要 标志 。 由 于 critical slowing down 效应 ( 详 
见 $1) ,临界 涨 落 在 实验 条 件 下 ,必然 会 处 在 非 平衡 
态 状态 ,因此 需要 构建 描述 涨 落 演化 的 动力 学 方程 ， 
并 将 其 与 流体 力学 方程 自 洽 地 耦合 起 来 。 

5) 临 界 涨 沙 的 粒子 化 :将 序 参量 ,或 者 重子 数 密 
度 的 涨 落 ,定量 转化 成 实验 可 观察 的 强 子 数 涨 落 , 需 
要 构建 合理 的 粒子 化 模型 。 

6) 数 据 分 析 : 贝 叶 斯 全 局 分 析 在 重 离子 碰撞 实 
验 中 有 着 广泛 应 用 。 对 于 寻找 QCD 临界 点 和 提取 
QCD 物质 有 限 密度 下 的 输 运 性 质 , 同 样 需要 拓展 已 
有 的 贝 叶 斯 分 析 工 具 。 


1 “理论 上 的 进展 


1.1 格 点 QCD 


由 于 符号 问题 , 格 点 QCD 不 能 直接 确定 QCD 
临界 点 是 否 存在 ,或 者 其 在 相 图 上 的 位 置 。 但 是 ,由 
于 格 点 可 以 第 一 性 原理 计算 重子 密度 较 低 区 域 (其 
体 而 言 ,ws/T < 2) 的 相 变 相关 属性 ,其 结果 是 寻找 临 
界 点 的 基石 。 

如 果 我 们 生活 在 一 个 夸克 质量 为 零 的 世界 AB 
么 手 征 对 称 性 是 精确 的 ,因此 ,重子 密度 为 零 时 ,从 
强 子 气体 到 QGP 的 相 变 伴随 着 手 征 对 称 性 从 破 缺 
到 还 原 的 二 级 相 变 。 在 真实 世界 ,夸克 质量 不 为 零 ， 
手 征 对 称 性 只 是 一 个 近似 。 相 应 地 , 相 变 是 没有 鲜 
明 分 界 的 平滑 相 变 , 相 变 温度 因而 无 法 严格 定义 。 


| 


文献 L5-7]。 此 外 , 重 离子 对 撞 物 理 近 一 两 年 部 分 有 


个 自然 的 选择 是 考虑 手 征 凝 聚 Cchiral 
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condesation ) 的 涨 落 ( 由 手 征 凝聚 susceptibility 刻画 ) 
随 温 度 的 变化 ,将 涨 落 峰值 对 应 的 温度 定义 为 相 变 
温度 。 通 过 泰勒 展开 , 格 点 不 仅 可 以 计算 pn = 0, 也 
可 以 计算 ws/ 不 太 大 条 件 下 的 手 征 凝 聚 ,进而 提取 
相 变 温度 。 由 此 得 到 的 相 变 温度 随 着 重子 化 学 势 的 
变化 ,可 以 写成 如 下 的 形式 : 


T.( He) Hp \2 HB 、4 
A nE ss 本 十 全 
E RAAE O 


式 中 :7 表示 重子 化 学 势 s = 0) 情 形 下 的 相 变 温 
度 ; 描 述 了 相 边 界 的 曲率 ,为 正 , 说 明 在 化 学 势 
较 小 的 时 候 , 相 变温 度 随 着 重子 密度 增加 而 降低 ;k， 
步 刻 画 有 限 密度 对 相 变 温度 的 修正 。HotQCD 
合作 组 "0 和 WB 合作 组 "得 到 的 结果 分 别 是 7, = 
(156.5 + 1.5) MeV, k, = 0.012(4), k, = 0.000(4) 和 
T, = (158.0 + 6) MeV ,k, = 0.015 3(18), k, =0.000 32 
(67), 代 表 了 这 方面 研究 的 最 新 进展 。 需 要 注意 的 
是 ,曲率 在 数值 上 较 小 ,四 阶 修正 几乎 可 以 忽略 ， 
说 明 相 边界 在 小 密度 区 域 可 以 近似 地 看 作 一 条 平行 
于 化 学 势 轴 的 曲线 , 见 图 2。 
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2 格 点 QCD 计算 得 到 的 手 征 相 变 温度 随 重子 化 学 势 
的 变化 5 


Fig.2 Chiral phase transition temperature vs. chemical 
potential from Lattice QCD" 


除了 手 征 凝聚 随 温度 变化 的 峰值 在 相 图 上 的 位 
置 , 格 点 还 可 以 算出 该 峰 对 应 的 宽度 。 在 临界 点 附 
近 ,宽度 会 变 窗 ,因为 对 应 的 峰会 变 陡 峭 。 文 献 L2] 
研究 发 现 ,一 直到 化 学 势 300 MeV ,宽度 几乎 不 随 化 
学 势 变 化 。 这 个 结果 意味 着 在 1 < 300 MeV 的 范 
内 ,临界 点 存在 的 可 能 性 非常 小 。 

除了 手 征 凝 聚 的 涨 落 , 格 点 还 可 以 系统 研究 守 
恒 荷 的 涨 落 ,比如 重子 荷 的 涨 落 。 这 些 涨 落 既 包括 


高 斯 涨 落 ,也 包括 非 高 斯 的 高 阶 涨 落 。 在 临界 点 附 
近 , 重 子 荷 的 涨 落 趋 于 发 散 , 且 非 高 斯 涨 落 发 散 得 更 
快 。 现 有 格 点 对 重子 荷 涨 落 的 计算 ?也 进一步 表明 
Ls < 300 MeV 的 范围 ,临界 点 并 不 存在 。 尤 其 重要 
的 是 ,在 这 个 范围 内 ,重子 荷 四 阶 涨 落 和 二 阶 涨 落 的 
比值 随 化 学 势 变化 的 格 点 结果 和 STAR 合作 组 的 实 
验 结果 符合 得 很 好 。 虽 然 这 个 比较 忽略 了 很 多 动力 
学 的 因素 ,但 它 表 明了 格 点 研究 对 实验 研究 相 变 可 
以 起 到 指引 作用 。 

从 更 广义 说 ,热力 学 量 的 涨 落 由 状态 方程 决定 。 
格 点 提取 有 限 密 度 状态 方程 的 主要 手段 是 基于 在 零 
密度 附近 的 泰勒 计算 展开 。 格 点 上 获取 有 限 密度 状 
态 的 进展 ,可 以 参考 文献 [12]。 
1.2 ”临界 涨 落 的 演化 

临界 点 附近 序 参 量 (the order parameter field) K 
波长 (大 尺度 ) 下 涨 落 会 被 增强 ,并 在 接近 临界 点 时 
随 着 关联 长 度 ( 的 发 散 而 趋 于 发 散 。 临 界 增强 的 
长 波长 涨 落 正 是 临界 乳 光 现象 背后 的 物理 机 理 。 具 
As Bl] QCD 临界 点 , 序 参量 是 手 征 凝 聚 和 能 量 ,重子 
密度 的 线性 组 合 。 因 此 , 序 参 量 的 涨 落 也 会 引发 末 
态 强 子 数目 Cmultiplicity) ,比如 质子 数目 的 涨 落 。 
基于 临界 点 附近 的 热力 学 普 适 行为 Cuniversality , 
非 高 斯 涨 落 , 比 如 三 阶 矩 Cskewness)、 四 阶 和 拢 
(kurtosis) 对 序 参量 关联 长 度 的 增长 更 加 敏感 。 理 
论 还 预言 ,在 临界 点 附近 ,四 阶 和 矩 随 化 学 势 将 展现 特 
定 的 非 单调 行为 。 该 非 单调 行为 被 认为 是 临界 点 存 
在 的 标志 性 信号 后 。 此 外 ,基于 并 合 模型 ,临界 涨 落 
也 可 能 通过 测量 轻 核 (light nuclei) APARIR IU, 

上 述 讨 论 基于 热平衡 条 件 下 的 临界 属性 。 涨 落 
和 耗 散 定理 (fluctuation- dissipation theorem) 告 诉 我 
们 , 涨 落 的 弛 驶 时 间 正 比 于 涨 落 的 幅度 。 更 具体 地 
说 ,临界 涨 落 的 弛 豫 时 间 m~e:。 这 里 临界 动力 学 指 
Bz > OOF QCD 临界 点 z ~ 3) ,表明 临界 弛 豫 时 
间 趋 近 于 无 穷 。 这 种 在 临界 点 附近 热平衡 过 程 变 慢 
的 现象 就 叫 critical slowing down 效应 。 由 于 火球 脱 
胀 冷却 和 critical slowing down 效应 ,在 重 离子 对 撞 
实验 中 , 序 参量 的 涨 落 最 终 在 劫难 逃 地 进入 非 平 衡 
状态 。 建 立定 量 描述 临界 涨 落 演化 的 理论 框架 是 寻 
找 临 界 点 需要 面 对 的 一 个 重大 挑战 。 

在 BEST collabration 成 立 之 前 ,对 于 非 平衡 临 
界 涨 落 已 经 有 了 一 系列 定性 研究 。 比 如 文献 [17] 就 
发 现 , 非 平衡 态 和 平衡 态 相 比 , 涨 落 不 仅 会 有 定量 
别 ,甚至 在 一 些 和 观测 量 紧密 相关 的 物理 量 上 ,有 定 
性 上 的 区 别 。 比 如 三 阶 矩 在 非 平 衡 态 下 的 符号 可 以 
和 平衡 态 下 的 期 待 值 相反 。 

QCD 临界 涨 落 演化 的 主要 定性 特征 可 以 通过 
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非 平衡 态 动 力学 中 的 Kibble-Zurek 机 制 来 描述 。 如 


的 数值 统计 量 。 但 是 其 挑战 在 于 描述 高 阶 涨 落 时 ， 


图 3 所 示 ,考虑 由 重 离子 对 撞 实 验 产生 的 火球 ,其 演 
化 对 于 相 图 上 一 条 经 过 临界 点 附近 的 轨迹 。 随 着 火 
球 靠近 临界 点 ,平衡 态 不 断 增 长 ,与 此 同时 ,达到 平 
衡 态 所 需 的 时 间 也 在 增长 ,因此 ,系统 实际 的 涨 落 最 
终 将 无 法 跟 上 平衡 态 的 增长 ,停滞 在 相 图 上 某 个 靠 
近 临 界 点 的 位 置 ,比如 图 中 的 O 点 。 因 此 ,系统 的 实 
际 涨 落 ,不 是 平衡 态 的 值 ,而 是 0O 点 附近 的 值 ,反映 
了 系统 对 O 点 的 记忆 。 这 个 效应 有 时 候 也 被 称 作 记 
忆 效 应 (memory effect)。 


Hy 


图 3 重 离子 碰撞 实验 产生 的 火球 在 临界 点 (红色 点 ) 附 近 的 
|2) 不 局 


LN 
Fig.3 Illustration of the non-equilibrium effect near the QCD 
critical point” 

BEST 合作 组 的 研究 极 大 推动 了 对 临界 涨 落 的 
定量 描述 。 从 理论 上 说 ,描述 非 平衡 态 临 界 涨 落 的 
基本 框架 是 涨 落 流体 力学 。 通 常 的 流体 力学 把 一 个 
多 体系 统 分 成 很 多 子 系统 ( 子 系统 有 时 也 被 称 作 流 
体 cell) ,描述 每 个 cell 的 平均 能 动量 ,以 及 守恒 荷 的 
时 间 演 化 ,空间 分 布 。 由 于 每 个 流体 cell 都 可 以 和 
周围 的 cell 发 生 能 动量 和 粒子 交换 ,因此 ,可 以 看 作 
巨 正则 系 综 发 生 涨 落 。 涨 落 流体 力学 的 目标 就 
是 在 描述 平均 能 动量 以 及 守恒 荷 变 化 之 外 ,进一步 
描述 其 涨 落 的 演化 。 如 果 对 涨 落 流 体力 学 输入 临界 
点 的 相关 信息 ,包括 临界 点 的 状态 方程 ,就 可 以 描述 
临界 "点 附近 的 动力 学 演化 。 

传统 上 描述 涨 落 流 体力 学 的 方法 是 “随机 方法 ” 
(stochastic approach) , EH Landau-Lifshitz 提出 的 办 
法 。 其 要 点 是 在 原 有 流体 力学 方程 中 加 入 “噪音 
项 ” ,噪音 的 幅度 大 小 可 以 通过 涨 落 - 耗 散 定理 与 输 
运 系 数 的 大 小 联系 起 来 。 最 近 几 年 ,一 个 新 的 方法 
被 引入 到 相对 论 涨 落 流体 力学 中 ,本 文 将 其 称 为 “ 决 
定论 方法 ”(deterministic appraoch)"”"。 其 主要 思想 
是 引入 一 系列 描述 涨 落 的 关联 函数 ,它们 和 通常 的 
流体 力学 自由 度 耦 合 在 一 起 ,满足 特定 的 


其 对 应 的 方程 形式 较 复 杂 。 

下 面 将 介绍 HYDRO+ 理 论 中 外, 它 是 “决定 论 方 
法 ”在 临界 点 附近 的 具体 应 用 。 为 了 描述 序 参量 MM 
的 长 波长 涨 落 6M,HYDRO+ 考 虑 了 如 下 的 Wigner 
函数 : 
Polt È) = | 7? < M (1,3 - p/2)0M (1,8 + 9/2) > 


(2) 
TA 5 (t,x) REE A Ex E ,高 斯 涨 落 在 尺 
度 1/0 下 的 强度 ,是 HYDRO+ 的 重要 自由 度 。 在 
HYDRO+ 中 ,go(t,x) 满 足 一 个 特定 的 驰援 方程 ,该 
方程 和 标准 的 能 动量 守恒 方程 和 重子 荷 守 恒 方 程 
起 ,构成 了 HYDRO+ 方 程 。 

需要 强调 的 是 , 涨 落 和 流体 力学 自由 度 是 耦合 
在 一 起 的 。 其 耦合 方式 主要 体现 在 :压强 不 再 只 是 
能 量 密度 和 守恒 荷 密度 的 函数 , 它 还 依赖 于 涨 落 
9o(t,*)。 这 是 因为 热力 学 函数 ,如 压强 ,在 临界 点 
附近 具有 临界 度 规 特征 (Critical scaling)。 这 些 特征 
源 自序 参量 的 涨 落 。 因 此 , 涨 落 的 非 平衡 态 演化 也 
会 反作用 于 热力 学 函数 ,所 以 在 HYDRO+ 中 ,压强 
也 是 g(t,x*) 的 泛 函 ,其 具体 形式 可 以 通过 统计 场 论 
的 方法 系统 推导 , 详 见 文献 [6]。 此 外 ,文献 [20-21] 
报道 了 将 Hydro+ 框 架 拓展 到 可 以 描述 非 高 斯 涨 落 
的 相关 进展 。 


2 ”模型 


2.1 ”初始 条 件 


如 前 所 述 , 和 LHC 以 及 RHIC 最 高 对 撞 能 量 相 
比 ,BES 能 量 下 火球 初始 状态 会 呈现 一 系列 新 的 特 
征 。 对 于 前 者 ,两 个 重 离子 核 相互 穿 过 的 时 间 
(Passing time) 是 0.15 fm 或 者 更 短 , 而 且 , 对 于 中 间 
R K OCRE y| < 1) ,近似 满足 longitudinal boost 
invariance。 这 些 特征 在 BES 能 量 下 并 不 成 立 。 对 
于 该 能 区 的 对 撞 , 穿 过 时 间 可 以 长 达 几 个 费 米 。 如 
能 量 密度 等 物理 量 对 快 度 更 加 敏感 ,使 得 
longitudinal boost invariance 不 再 是 一 个 很 好 的 近 
似 。 由 此 可 见 ,BES 能 区 的 初始 条 件 模型 构建 复杂 
度 较 高 。 

BEST 合 作 组 成 员 在 初始 条 件 模型 方面 一 个 重 
要 工作 是 发 展 了 “动态 初始 化 ”方法 ”。 如 前 所 述 ， 
TE BES 能 区 ,对 撞 离 子 穿 过 时 间 可 以 长 达 几 个 费 


“deterministic” 方 程 。 研 究 表 明 , 该 方法 成 功 描述 
系列 重要 的 非 平衡 态 效 应 。 该 方法 的 优点 在 于 对 计 
算 资 源 的 要 求 远 低 于 随机 方法 ,因为 后 者 需要 极 高 


Ko 意味 着 在 对 撞 的 两 个 核 在 相互 穿 过 之 前 ,火球 
中 已 经 有 一 部 分 趋 于 热 化 ,可 以 被 流体 力学 描述 。 
“动态 初始 化 ?方法 将 没有 热 化 部 分 的 贡献 看 作 源 


040010-4 


202306.00136v1 


chinaXiv 


核 技 术 


ChinaXiv 合 作 期 刊 


2023, 46: 040010 


(source) ,对 已 经 热 化 部 分 求解 含 源 的 流体 力学 方 
程 。 文 献 [9] 将 该 方法 得 到 的 初始 条 件 输入 到 流体 
力学 的 模拟 中 ,并 与 已 有 的 数据 进行 了 比较 ,包括 质 
子 数 对 快 度 的 依赖 。 这 些 研究 验证 了 “动态 初始 化 ” 
方法 是 一 个 在 唯 象 上 较为 成 功 的 模型 。 
于 初始 状态 重子 密度 空间 分 布 的 演化 会 对 末 
态 质 子 涨 落 产生 贡献 。 从 寻找 临界 点 的 角度 看 ,这 
些 贡献 与 临界 点 无 关 , 是 需要 扣除 的 本 底 。 因 此 , 现 
有 初始 化 模型 需要 进一步 完善 对 初 态 涨 落 的 定量 
描述 。 
2.2 ”包含 了 临界 点 的 状态 方程 

状态 方程 是 流体 力学 的 重要 输入 ,因此 ,为 了 求 
解 临界 点 附近 的 流体 力学 方程 ,首先 需要 构建 包含 
临界 点 的 状态 方程 。QCD 临界 点 被 认为 和 3D Ising 
模型 属于 同一 个 universality class。 因 此 ,构建 状态 
方程 需要 满足 两 个 重要 约束 :1) 在 临界 点 附近 和 3D 
Ising 模 型 具有 可 以 一 一 对 应 的 临界 属性 ;2) 在 小 重 
子 密度 区 和 格 点 的 结果 一 致 。 此 外 ,状态 方程 在 其 
适用 的 相 图 范围 ,需要 满足 热力 学 稳定 性 和 因果 律 
(声速 低 于 光速 ) 等 限制 。 显 然 ,即便 加 上 这 些 限 秆 
构建 状态 方程 的 方法 也 不 唯一 。BEST 合作 组 构造 
的 状态 方程 发 表 于 文献 [23j 中 。 由 于 QCD 临界 点 
和 3D Ising 模型 之 间 的 映射 依赖 于 几 个 非 普 适 性 的 
参数 ,因此 ,BEST 状态 方程 也 依赖 于 这 些 参数 。 具 
体 而 言 ,它们 包括 临界 点 的 位 置 (7.,y),Ising 模 型 
相 图 上 的 温度 轴 和 磁场 轴 映 射 到 QCD 相 图 之 后 和 
QCD 相 图 的 温度 ,化 学 势 轴 的 夹 角 ai、a, 以 及 参数 
(w,p) 控 制 临 界 区 域 在 温度 、 密 度 方向 的 范围 ,这 些 
参数 刻画 了 QCD 临界 点 的 特征 。 通 过 实验 和 以 
BEST 状态 方程 为 输入 的 动力 学 模拟 计算 比较 ,有 可 
能 限制 甚至 确定 这 些 参数 。 

对 于 寻找 临界 点 ,BEST 状态 方程 的 应 用 较为 ) 
泛 。 首 先 , 它 可 以 作为 求解 流体 力学 的 输入 ;其 次 ， 
它 可 以 用 来 研究 平衡 态 条 件 下 ,临界 涨 落 对 于 重子 
化 学 势 的 依赖 。 比 如 文献 [24] 分 析 了 重子 密度 四 阶 
涨 落 ,发 现 如 果 只 考虑 临界 普 世 性 的 效应 ,四 阶 涨 落 
会 随 着 化 学 势 增长 先 下 降 , 再 上 升 到 峰值 。 但 是 如 
果 考 虑 了 非 普 世 的 贡献 ,下 降 的 行为 有 可 能 会 消失 ， 
只 有 上 升 到 峰值 的 表现 不 会 变化 。 

需要 补充 的 是 ,求解 流体 力学 还 需要 输入 输 运 
系数 ,包括 帮 滞 系数 、 弥 散 常 数 。 这 些 系数 很 难 通 过 
格 点 计算 提取 ,此 外 输 运 系数 也 有 特定 的 临界 表现 。 
参数 化 输 运 系数 在 BES 能 区 的 表现 ,并 包含 正确 的 
临界 行为 ,这 是 一 个 需要 完善 的 工作 。 


Ls 
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2.3 ”流体 力学 演化 

BEST 合 作 组 发 展开 发 了 两 套 独立 的 三 维 流体 
力学 求解 代码 :1) 将 麦 吉尔 大 学 研究 组 原 有 的 
MUSIC 代码 拓展 到 了 BES 能 区 (以 下 简称 
MUSIC) ; 2) 俄 辫 俄 州立 大 学 研究 组 开发 的 
BESHydro 代码 31( 在 BEST 之 外 ,CCNU 的 研究 组 
最 近 开 发 了 能 在 GPU 上 并 行 的 代码 CLVIS29) 。 文 
献 L12] 报 道 了 这 两 个 不 同 代 码 之 间 的 相互 验证 。 

在 输入 了 初始 条 件 、 状 态 方程 和 输 运 系数 之 后 ， 
流体 力学 可 以 给 出 对 撞 产 生 的 火球 在 相 图 的 轨迹 ， 
定量 地 回答 了 重 离子 对 撞 实 验 扫描 了 相 图 的 哪些 区 
域 。 图 4 展示 了 19.6 GeV 下 三 个 有 代表 性 的 空间 快 
度 所 对 应 的 轨迹 。 值 得 注意 的 是 , 当 快 度 从 中 心 区 
域 jn,| < 0.5 变化 到 较 大 值 时 0.5 < y, < 1.5, 对 应 的 
重子 密度 变化 率 50~100 MeV。 这 意味 着 可 以 把 扫 
苗 快 度 作 为 扫描 相 图 的 补充 手段 。 
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—— With Critical Point 
0.25 LL-===: No Critical Point 
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Do 
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0.05 F 
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“Its GeV 


4 通过 求解 流体 力学 得 到 的 火球 在 不 同 中 心 快 度 办 在 相 

上 所 对 应 的 轨迹 站 

Fig.4 Trajectories on the QCD phase diagram from 
hydrodynamic simulations for BESII" 


基于 包括 了 临界 点 的 状态 方程 和 流体 力学 , 文 
献 [1 进 一 步 分 机 了 临界 点 对 火球 轨迹 的 影响 。 图 
4 中 实 线 、 虚 线 分 别 表示 状态 方程 中 包含 了 和 没有 
包含 临界 点 时 对 应 的 轨迹 。 从 中 不 难 发 现 , 临界 点 
对 火球 轨迹 的 影响 并 不 显著 。 进 一 步 计算 还 发 现 ， 
末 态 粒子 的 分 布 也 没有 受到 临界 点 的 显著 影响 。 这 
说 明 涨 落 观 察 量 ,而 不 是 粒子 分 布 ,对 临界 点 更 
敏感 。 

对 于 有 限 密 度 的 QCD 物质 ,重子 密度 演化 是 流 
体力 学 演化 的 重要 组 成 部 分 。 重 子 密 度 弥 散 这 一 重 
要 的 物理 过 程 被 输 运 系数 重子 弥散 常数 控制 。 文 献 
[27] 研 究 了 末 态 质子 数 分 布 对 重子 弥散 常数 的 依赖 
关系 。 从 物理 上 说 ,弥散 过 程 会 使 末 态 分 布 和 初始 
重子 数 分 布 不 同 ,因此 ,有 可 能 通过 实验 提取 弥散 常 
数 。 如 图 5 所 示 , 质 子 数 在 快 度 空 间 的 分 布 确实 对 
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弥散 常数 有 一 定 的 敏感 性 , 实 线 和 虚线 分 别 对 应 重 F, sea een aes gs be pa 

子 弥散 常数 为 零 和 不 为 零 的 结果 。 但 是 ,由 于 现在 献 研究 了 其 在 唯 象 上 对 观测 量 的 影响 ,结论 是 :由 于 
对 初 态 重子 数 分 布 的 认 知 有 限 ,因此 ,在 现 阶段 , 直 重 离 子 对 撞 实 验 的 多 重复 杂 因素 ,观测 量 对 弥散 党 
接 基于 流体 力学 模型 提取 该 弥散 常数 还 有 很 大 的 不 。”” 数 的 临界 行为 敏感 度 非常 有 限 。 
确定 度 ,但 这 是 一 个 值得 努力 的 方向 。 在 临界 点 附 
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Fig.5 Proton multiplicity as a function of momentum rapidity™ 


2.4 ”临界 涨 落 的 演化 

§2.3 讨论 了 没有 考虑 涨 落 的 流体 力学 演化 ,本 
临界 涨 落 演化 的 进展 。 
基于 Hydro+ 框 架 ,文献 [29-30] 研 究 了 在 简化 
的 类 似 重 离子 对 撞 实验 的 条 件 下 ,考察 了 非 平衡 态 
高 斯 涨 落 go (t,x*)( 定 义 见 $2) 的 演化 ,获得 了 如 下 有 
价值 的 信息 。 
首先 ,这 些 工作 验证 了 平衡 态 涨 落 和 非 平衡 态 


节 介 


ASS 


涨 落 的 演化 是 不 一 样 的 。 图 6 中 ,不 同 颜色 的 曲线 i he eee 
对 应 不 同 的 时 间 点 ,平衡 态 和 非 乎 衡 态 涨 落 分 别 由 (a) 系统 靠近 临界 点 ,(b) 系统 远离 临界 点 

实 线 和 虚线 画 出 。 如 图 6 所 示 , 在 系统 开始 靠近 临 Fig.6 Non-equilibrium a the critical non-Gaussian 
界 点 的 时 候 , 非 平衡 态 涨 落 ( 实 线 ) 努 力 追 赶 不 断 增 (a) Near the critical point, (b) Away from the critical point 
长 的 平衡 态 涨 落 ( 虚 线 ) ,但 被 critical slowing down zx 帝 度 。 虽 然 序 参量 涨 落 这 一 个 自由 度 的 贡献 在 临界 


效应 压低 。 当 系统 开始 远离 临界 点 时 ,平衡 态 涨 落 。 点 附近 会 被 加 强 ,但 是 QCD 物质 还 存在 大 量 和 胶 
会 衰减 ,但 是 非 平 衡 态 涨 落 衰减 较 慢 ,能 够 保存 部 分 子 ,夸克 相关 的 自由 度 ,后 者 压低 了 序 参量 涨 落 的 反 


关于 临界 友 的 “记忆 ,参看 前 文 。 这 也 说 明 , 我 们 需 。 作用。 因此 对 于 寻找 QCD 临界 点 的 模型 计算 ,可 以 
要 通过 分 析 非 平衡 态 涨 落 ,提取 临 春 信和 号。 先 忽略 这 些 反作用 。 
其 次 ,守恒 律 对 非 平衡 态 演化 有 非常 重要 的 作 如 前 所 述 ,， 基于 “deterministic 方法 ”的 


用 。 由 于 重 了 密度 是 守恒 的 ， 因此 长 波长 (小 O) 的 se BEST f 
ee E. ENHET BERROETA TE T, SMDS 
FER PASEAR REAREA 通过 求解 包含 临界 点 的 随机 流体 力学 方程 ,计算 了 


的 影响 可 以 忽略 。 这 里 反作用 是 指 在 涨 落 处 在 非 疗 ”高 斯 和 非 高 斯 涨 落 的 演化 ,结果 表明 :该 方法 也 可 以 


in 


衡 态 条 件 下 ,流体 力学 自由 度 ,如 能 量 密度 ,流速 度 ， ee et ne . 
它们 的 演化 会 发 生 怎样 的 改变 。 文 献 [7-8] 的 结果 We oe A 
发 现 , 涨 落 演化 对 能 量 密度 ,流速 度 的 影响 最 大 不 超 25 ”粒子 化 (particlization) 

过 0.1% 这 个 量 级 。 对 于 这 个 结果 ,一 个 直观 的 解释 将 流体 力学 的 结果 和 实验 观察 量 联系 起 来 的 关 
是 : 广 延性 热力 学 函数 的 大 小 正比 于 系统 自由 度 的 键 步 又 是 粒子 化 (particlization)。 它 的 目标 是 将 能 
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量 密度 ,重子 密度 以 及 流速 度 (fluid velocity) 的 时 空 
分 布 ,以 及 它们 的 涨 落 ,和 实验 可 以 测量 的 强 子 产 
额 ,动量 空间 的 分 布 以 及 相应 的 涨 落 , 定 量 地 关联 起 
来 。 因 此 ,对 临界 涨 落 进行 合理 的 粒子 化 ,是 寻找 临 
界 点 必须 解决 ,又 充满 挑战 性 的 问题 。 

从 信息 论 的 角度 而 言 , 粒 子 化 问题 的 难度 在 于 
粒子 动量 空间 分 布 函数 包含 了 比 流体 力学 解 更 多 的 
信息 。 更 具体 地 说 ,对 分 布 函数 jj 进行 积分 可 以 得 到 
能 动量 以 及 守恒 和 荷 ,因此 ,后 者 是 对 分 布 函数 是 一 种 
限制 。 但 是 有 无 穷 多 的 可 以 满足 这 些 限 制 ,因此 ， 
唯一 确定 分 布 函 数 ,需要 额外 地 输入 或 者 假设 。 在 
不 考虑 涨 落 的 条 件 ,被 广泛 使 用 的 Copper-Frye # 
化 方案 假设 了 分 布 函数 是 热平衡 态 分 布 。 那 么 通过 
能 量 、 守 恒 蓓 确定 相应 的 温度 、 化 学 势 , 对 应 的 分 布 
函数 也 就 可 以 确定 下 来 。 但 是 这 个 方法 显然 不 能 
接 推广 到 了 粒子 化 处 在 非 平 衡 状 态 下 的 临界 涨 落 。 

文献 [32] 假 设 强 子 分 布 的 涨 落 是 源 于 序 参量 场 
的 涨 落 , 后 者 由 于 和 强 子 场 耦合 ,引发 强 子 质量 的 涨 
落 , 进 而 带 来 粒子 数 分 布 的 涨 落 。 在 这 样 的 假设 下 ， 
文献 [19] 将 Hydro+ 模 拟 的 结果 转化 成 粒子 数 (如 质 
子 数 ) 的 高 斯 涨 落 。 这 个 方法 也 可 以 相对 直接 地 推 
广 到 高 阶 ( 非 高 斯 ) 涨 落 。 其 缺点 在 于 对 于 粒子 数 涨 
落 和 序 参 量 涨 落 的 微观 原因 ,有 一 个 较 强 的 假设 ,该 
假设 还 没有 办 法 被 验证 。 与 此 平行 ,文献 L33] 拓 展 
了 “局 域 守恒 ”假设 ,发 展 了 不 同 的 粒子 化 涨 落 方法 。 
最 近 ,文献 L34] 基 于 信息 论 ,提出 了 新 的 解决 方法 ， 


有 一 定 启发 性 。 粒 子 化 临界 涨 落 ,依然 是 一 个 开放 
性 的 问题 。 


2.6 ”其 他 进展 


这 里 将 简 述 BEST 合 作 组 一 些 其 他 相关 进展 。 

临界 点 以 及 一 级 相 变 点 存在 不 仅 会 影响 火球 在 
流体 力学 阶段 ,也 会 影响 其 在 强 子 气体 态 阶段 的 演 
化 。 此 前 描述 后 者 的 输 运 模型 并 没有 考虑 临界 效 
应 。 文 献 的 作者 通过 引入 平均 场 , 发 展 了 高 密 区 核 
物质 状态 方程 模型 (Vector Density Functional, VDF) 
模型 。 该 模型 的 优点 是 可 以 直接 和 强 子 输 运 的 模拟 
代码 SMASH 融合。 文献 [15] 的 研究 还 表明 aR 
气体 阶段 的 临界 效应 也 潜在 对 观测 量 有 影响 。 关 于 
输 运 模型 的 近期 进展 ,可 以 参考 文献 [36j]。 

此 外 ,BEST 合 作成 员 还 以 提取 有 限 密度 强 相互 
作用 物质 密度 和 临界 点 信息 为 目标 ,开发 相应 的 全 
局 贝 叶 斯 分 析 程 序 。 由 于 模型 参数 的 增加 ,实现 这 
些 目标 的 所 需 计 算 资源 比 此 前 对 高 能 的 重 离子 碰撞 
试验 结果 进行 分 析 所 需 的 资源 可 能 高 一 到 两 个 数 
量 级 。 


3 ”总 结 与 展望 


从 重 离子 对 撞 实验 中 提取 有 限 密度 区 强 相互 物 

质 的 属性 ,包括 寻找 其 可 能 存在 的 相 边 界 和 临界 点 ， 

需要 构建 一 个 复杂 的 ,可 以 描述 火球 演化 的 动力 学 
框架 。BEST 合 作 组 的 工作 ,对 这 个 框架 中 最 重要 的 
几 个 部 件 的 构建 ,都 起 到 不 可 或 缺 的 贡献 。 它 们 
包括 : 

1) 一 个 复杂 但 可 以 在 数值 上 实现 的 ,包含 了 重 
子 stopping 和 有 限 passing time 的 初始 状态 模型 。 

2) 将 此 前 已 有 的 相对 论 流体 力学 模型 推广 到 
BES 对 接 能 量 ,并 可 以 用 于 模拟 包括 重子 荷 ,奇异 荷 
在 内 的 守恒 荷 演 化 。 

3) 构 建 了 包含 
学 模型 的 输入 。 

4) 发 展 了 可 以 演化 临界 涨 落 的 相对 论 涨 落 流体 
力学 理论 ,包括 HYDRO+ 理 论 , 并 初步 应 用 于 数值 
模拟 中 。 

5) 探 索 了 粒子 化 临界 涨 落 的 方法 ,为 定量 描述 
实验 观察 量 打下 基础 。 

上 述 这 些 “ 部 件 ”, 已 经 被 世界 各 地 的 研究 者 广 
泛 应 用 到 了 相关 研究 中 去 。 对 于 RHIC 能 量 扫描 相 
关 物 理 的 研究 ,正在 进入 从 定性 到 定量 的 阶段 。 

展望 未 来 ,世界 的 主要 科技 强国 ,包括 德国 、 俄 
罗斯 和 中 国 , 正 在 计划 或 筹备 在 大 型 科学 装置 上 的 
重 离子 对 撞 实 验 。 这 些 实验 的 主要 目标 就 是 进一步 
研究 中 高 密度 区 的 QCD 相 图 。 其 中 低温 高 密 核 物 
质谱 仪 (CSR-External-target Experiment, CEE) 于 
2019 年 获 基金 委 国 家 重大 科研 仪器 研制 项 目 资助 。 
CEE 将 是 我 国 第 一 台 运 行 于 GeV 能 区 的 、 自 主 研制 
的 、 基 于 HIRFL-CSR 和 HIAF "的 大 型 核 物理 实验 
AB. CEE 实验 对 寻找 一 级 相 变 的 信号 ,以 及 确定 
QCD 临界 点 在 特定 重子 密度 区 是 否 存在 ,有 着 极其 
重要 的 意义 。 这 些 实验 ,按照 设计 ,可 以 探测 比 
RHIC-BES 实验 重子 密度 更 高 的 区 域 , 也 将 对 动力 
学 模型 的 构建 提出 新 的 要 求 和 挑战 。 一 方面 ,BEST 
合作 组 取得 的 结果 ,可 以 作为 未 来 这 些 研究 的 起 点 ; 
另 一 方面 ,一 些 BEST 合 作 组 没有 解决 的 难题 ,需要 
在 未 来 被 解决 ,它们 包括 :如 何 描述 一 级 相 变 区 域 的 
流体 力学 以 及 相关 涨 落 的 演化 ;如 何 将 现 有 框架 的 
适用 范围 推广 到 更 低 的 对 撞 能 量 ( 比 如 5 GeV) 以 
下 。 对 于 这 些 碰撞 ,处 在 流体 力学 可 以 适用 的 阶段 
相对 较 短 ,因此 呼唤 新 一 代 动 力学 模型 的 产生 和 
发 展 。 
作者 贡献 声明 本 文 唯 
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